We claim: 



1 . A method for inferring a network relationship between genes, comprising: 

(a) providing a quantitative time course data library for a set of genes of 
an organism, said library including expression results based on time course of 
expression of each gene in said set of genes, quantifying an average effect and 
measure of variability of each time point on each other of said genes; 

(b) creating a sparse matrix from said library, said matrix having zero 
coefficients removed therefrom; 

(c) generating a set of linear differential equations from said matrix; and 

(d) solving said set of equations to produce said network relationship. 

2. The method of claim 1, wherein said zero coefficients are identified using 
Akaike's Information Criterion (AIC). 

3. The method of claim 1, wherein said differential equation is 

x ( t) 

in which the vector — v } contains the amount of expressed cDNA as a function of 

A . , . ,-i 

time, and the matrix — is a constant with units second . 

4. The method of claim 1, wherein said matrix contains elements Ay 9 wherein 
Ay represents the effect of gene yon gene i, and wherein [a^] represents the 

reaction time for said effect of gene j on gene L 

5. The method of claim 1, wherein said differential equation solved is 

Attorney Docket No: GENN1009US1 DBB Express Mail No.: EV327622585US 

dbb/genn/1 009US 1 /GENN. 1 009US 1 .doc 



22 



2i(t) = exp 



v 



6. The method of claim 1, wherein said exponent At (exp(A)) is solved using the 
formula: 



exp A_ 



00 1 • 
) i = 0 1 ' 



7. The method of claim 1, wherein said differential equation is estimated by 
solving the difference equation: 

Ax 

8. The method of claim 1, wherein said sparse matrix further comprises an error 
estimated using the formula: 

x(r+Ar) -x(0 = A/*Avr(r) -f_e(r). 

9. The method of claim 8, wherein said error has a normal distribution 
independent of time according to the formula: 



/(iKOio 2 ) 



T 



W27IG 2 ; 



ex 



2a 2 



wherein standard deviation o is equal for each of said genes at all times. 

10. The method of claim 1, wherein the maximum likelihood estimate of the 
variance a 2 is determined by maximizing the log-likelihood function with respect to 
a 2 using the formula: 
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1 1 . The method of claim 1 0, wherein said variance a is determined using the 
formula: 



*2 1 ^ 
a = — > 

i = l 



Z t T 

( *T -^f- 1) * ( ^ - 1) + ( '/ - '/ - 1) d- - 1 



12. The method of claim 2, wherein said AIC is minimized using the formula: 

a m o T '°g- likeliho ° dofthe "I, 0 f number of estimated 1 n ~ 
AIC-2 |^ estmated model J +2 [ parameters J (16) 

M r A 

13. The method of claim 1, wherein mask — is used to set matrix elements of — 



equal to zero using the formula: 

A = M°A, 

M 

where o denotes an element-wise product, and mask — is a matrix whose elements 



are either one or zero. 

1 4. The method of claim 1 3, wherein matrix elements are set to zero by applying 
a mask M produced by minimizing the formula: 
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A 

thereby yielding the maximum likelihood estimate — 



15. The method of claim 2, wherein said AIC is minimized according to the 
formula: 



AIC= h m In 



2 no 



+ nm + 2 # ( 1 + [ sum of the maskelements M r ,yj) , 



M 

1 6. The method of claim 13, wherein said mask — is selected to minimize AIC 



calculated using the formula: 



AIC= n /win 



27ia 



«m + 2 # ( 1 + [ sum of the maskelements M^) , 



17. A medium containing one or more results of network relationships between 
genes calculated using a method of claim 1 stored thereon. 

18. A method for determining the statistical significance of network relationships, 
comprising: 

(a) calculating the average log-ratio of expression for each gene at each 
time point; 

(b) calculating the standard deviation from all measurements; 

(c) calculate the joint probability; and 

(d) adopting a criterion for statistical significance. 
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19. The method of claim 18, wherein said step (a) is determined using the 
formula: 

20. The method of claim 18, wherein step (b) is determined using the formula: 



A4 =<i± I (*/.W) 2 > 

1 / = I A' - 1,2 



Xj t [k] 



in which J " is the data value of measurement k at time point i for gene j. 

21. The method of claim 18, wherein the joint probability for j to be larger in 
absolute value than the measured values . is calculated using the formula: 



1 = 1 i = 1 











> 


x ji 



n 

/ = i 



1 - erf 



wherein ef f is an error fiinction. 

22. The method of claim 18, wherein a significance level a is selected. 

23. The method of claim 18, wherein the null hypothesis is rejected if Pi < a. 
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24. The method of claim 18, wherein the null hypothesis is rejected if P < a> 
wherein n is the number of time points at which gene expression is evaluated. 

25. A method for determining the statistical significance of network relationships, 
comprising: 

(a) calculating the average log-ratio of measurements of expression for 
each gene at each time point using the formula: 

(b) calculating the standard deviation of said measurements using the 
formula : 



a 



> / = 1 A = 1,2 



in which is the data value of measurement k at time point / for gene j. 

(c) calculating a joint probability for % to be larger in absolute value 

than measured values jj i calculated using the formula: 



'■n'( = n'< 











> 





= n 

i = i 



1 - erf 
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wherein er f is an error function; and 

(d) applying a criterion for statistical significance to determine whether a 
null hypothesis is rejected. 

5 26. The method of claim 25, wherein the null hypothesis is rejected if P < a, 
wherein n is the number of time points at which gene expression is evaluated. 
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